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Dynamical chaos has recently been shown to exist in the Gaussian approximation in quan- 
tum mechanics and in the self-consistent mean field approach to studying the dynamics of 
quantum fields. In this study we first show that any variational approximation to the dynam- 
ics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian 
dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear 
and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact 
quantum evolution. We then restrict attention to a system of two biquadratically coupled 
quantum oscillators and study two variational schemes, the leading order large N (four 
canonical variables) and Hartree (six canonical variables) approximations. The chaos seen 
in the approximate dynamics is an artifact of the approximations: this is demonstrated by 
the fact that its onset occurs on the same characteristic time scale as the breakdown of the 
approximations when compared to numerical solutions of the time-dependent Schrodinger 
equation. 

PACS numbers: 05.45. +b, 03.65. Sq, 2.30 Wd, 03.65 -w 



I. INTRODUCTION 

There are many situations in quantum mechanics and 
field theory where one hopes that one dynamical degree 
of freedom can be considered "classical" or "semiclas- 
sical." In the dynamics of the early universe, one usu- 
ally imagines that gravitational energy can be transferred 
to particle production, with the gravitational field being 
treated semiclassically, i.e., the quantum matter fields 
evolve in a background "classical" gravitational field, the 
dynamics of which is in turn determined from the ex- 
pectation value of the energy momentum tensor of the 
quantum field. Similarly in pair production from strong 
electric fields, one attempts to describe the background 
electric field "classically" and solve for the dynamics of 
the quantum degrees of freedom in this background field. 
The time dependence of the electric field is governed by a 
Maxwell equation in which the right hand side is the av- 
erage value of the current of the produced pairs. In this 
sort of approximation of a quantum system coupled with 
a semiclassical degree of freedom such as a coherent elec- 
tric or gravitational field, the approximate dynamics of 
the quantum system can become chaotic. This was first 
described by us, and termed "semiquantum chaos" [§],§||. 
A closely related result having the same cause is "semi- 
quantal chaos" || which occurs in the time-dependent 
Gaussian approximation for the dynamics of quantum 
systems. 

What happens in these dynamical approximations is 
that the time evolution of the parameters governing the 
shape of the quantum mechanical wave function (or func- 
tional) becomes sensitive to the initial conditions. In 
this paper we will first establish that this behavior can 



occur in any variational approximation to the quantum 
dynamics {e.g., time dependent Hartree approximation). 
We will then focus on exactly the same model system 
treated in Refs. |j]J|, namely a system of two coupled 
oscillators described by the Lagrangian: 

L = ii 2 + ii 2 -i ( ™ 2 + e 2 AV • (1) 

This system of two nonlinearly coupled oscillators arose 
from studying the problem of pair production of charged 
pions in a strong external electric field Q (quantum fluc- 
tuations of the electric field were ignored). In momen- 
tum space, the individual modes of the pion field dis- 
played chaotic behavior. The two-oscillator problem re- 
sults from ignoring all but the k = mode for the quan- 
tum field. In the Lagrangian ([!]), the A oscillator repre- 
sents the time-dependent electromagnetic field and the x 
oscillator, the k — mode of the charged pion field. 

Treating the electromagnetic (A) field classically is the 
standard first term in a large TV expansion D and is re- 
lated to the classic problem treated first by Schwinger 
Q on pair production from external fields. Because such 
semiclassical methods are often used in initial value prob- 
lems in quantum field theory, we hope to understand 
the origin of the chaos by considering a simple quantum 
mechanical model. To this must be added the impor- 
tant point that while accurate numerical solutions to the 
quantum mechanical problem are available to test the va- 
lidity of approximations, such a luxury is not available in 
field theory. 

The semiclassical calculation is equivalent to a Gaus- 
sian variational approximation to the field theory (see 
Ref. for more details and an explanation of dissi- 
pation and decoherence in this approximation). As we 
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show later, all variational approximations to quantum 
dynamics lead to classical Hamiltonian dynamics for the 
variational parameters (the Gaussian approximation is a 
special case of this general result). Therefore, since the 
resulting Hamiltonian dynamics is gencrically nonlinear, 
chaos can be present in the approximate dynamics. 

We will consider two variational approximations which 
are equivalent to two different assumptions about the 
fluctuations of the A oscillator. The first approximation 
(semiclassical or leading order large N) is the assump- 
tion that we can ignore all quantum fluctuations of the 
A oscillator (the quantum mechanical version of the elec- 
tromagnetic field). This is equivalent to assuming 

(A 2 *) = (A) 2 {x). (2) 

The Hartree approximation assumes that we include only 
Gaussian fluctuations of both quantum oscillators. This 
assumption implies a factorization of the expectation val- 
ues as: 

(A 2 x) = (A 2 )(x). (3) 

Both of these approximations, since they include only 
Gaussian fluctuations, are really no different than a par- 
ticular phase space ensemble of classical solutions of the 
equations of motion for the coupled classical oscillator 
problem, with a particular initial condition implementing 
the uncertainty relation. Consequently, the same chaos 
discussed above will also be found in the corresponding 
dynamics of the classical Liouville equation when only 
Gaussian fluctuations are allowed. This aspect of the 
Gaussian approximation we will discuss elsewhere Q . 

Our numerical results show that in the Hartree approx- 
imation, the onset of chaos, as a function of parameters 
of the Hamiltonian, is marginally delayed as compared to 
the large N (semiclassical) approximation. We find that 
both approximations diverge from the exact numerical 
simulation of the Schrodinger equation at approximately 
the same time. After that time, the Hartree approxima- 
tion qualitatively tracks the general features of the exact 
simulation better than the large N approximation. By 
direct comparison with the exact numerical solution we 
also find that chaos in the variational approximations 
occurs roughly on the same time scale as when these 
approximations diverge from the exact numerical solu- 
tion. As is known on general grounds 0], expectation 
values of the full quantum system (which are the vari- 
ational parameters of the classical Hamiltonian dynam- 
ics) are insensitive to initial conditions: Our results are 
completely consistent with this fact. Our interpretation 
of the above results is in accord with that of Sundaram 
and Milonni |ic| who have argued that the chaos seen 
in the approximate dynamics is not a fundamental fea- 
ture of the full quantum dynamics but simply reflects a 
breakdown of the approximation scheme, ft was further 
argued in Ref. ]Io[ that the approximations are unreliable 
when either the classical equations are already chaotic or 
when the approximate dynamics is chaotic. To test the 



second part of this statement we explored nonchaotic pa- 
rameter regimes for the approximate dynamics (but not 
too far from the onset of chaos) and found essentially no 
improvement in the agreement between the exact quan- 
tum and approximate calculations. Thus the existence of 
chaos is insufficient to assess the accuracy of the approx- 
imations: apparently the breakdown time (in terms of 
natural time scales) is the same whether chaos is present 
or not. 

It is strong nonlinearity rather than just chaotic dy- 
namics which leads to the breakdown of the approxima- 
tions. This is hardly surprising: the Gaussian approxi- 
mations are equivalent to truncating a cumulant expan- 
sion at second order. If the exact dynamics is strongly 
nonlinear, higher order cumulants are generated and a 
second order truncation quickly becomes invalid. For 
chaotic systems, Ref. [[it)) provides a simple analytic ar- 
gument, but the statement is true more generally. 

These results might seem to put very strong constraints 
on dynamical mean field approximations in quantum field 
theory, especially at strong coupling. However, the field 
theoretic analog to the above problem has a very large 
(formally infinite) number of degrees of freedom. For 
example, in the field theoretic case, Eqn. ([!]) is a radical 
truncation of the full Lagrangian 

L = \(d M - ieA^\ 2 - 1(^A, - d v A^f - m 2 ^<\> . (4) 

In the leading order large N approximation, A is still 
treated classically but it is now coupled to a very large 
number of fluctuating degrees of freedom. Before defini- 
tive statements regarding the accuracy of mean field ap- 
proximations can be made, two issues have to be clari- 
fied. The first has to do with the fact that even though 
the individual trajectories of the Fourier modes may 
be chaotic and far from the exact solution, what really 
matters is the summed contribution (i.e., the statistics 
of the distribution of trajectories) and this may have a 
much more benign character. The second point is related 
to the onset of chaos as the number of degrees of freedom 
is varied. The importance of this question was noted by 
Ford |y]] but it has not been studied in any detail in 
the literature. Thus, it is still an open question whether 
chaos in the mean field approximation in field theory is 
as serious an obstruction as suggested by the study of 
lower dimensional systems. 

The rest of the paper is set out as follows. First (Sec. 
II) we prove the general result that all variational approx- 
imations lead to a Hamiltonian dynamics for the varia- 
tional parameters. Then, in Sec. Ill wc explicitly discuss 
the Hamiltonian dynamics for the two oscillator prob- 
lem, both in the large N and Hartree approximations. 
In Sec. IV we briefly describe our numerical approach 
to the exact solution of the two coupled oscillator prob- 
lem. We then compare numerical simulations of the two 
variational approximations with the exact solution of the 
Schrodinger equation. Finally, in Sec. V, we state our 
conclusions and discuss the implications of our results. 
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II. THE TIME-DEPENDENT VARIATIONAL 
PRINCIPLE 

The Schrodinger equation can be reduced to a system 
of ordinary differential equations for some variational pa- 
rameters by constraining the wave function to be of a 
particular form. In this section we show that any vari- 
ational calculation of the wave function will lead to a 
Hamiltonian dynamics for the variational parameters. 

The starting point for a variational calculation is 
Dirac's action principle [jl2| which can also be used to 
derive the Schrodinger equation as shown below. We be- 
gin by defining the action: 



S 



dt (WiS- -i2]*)/(*|*) 



ti 



dt 



The time dependent Schrodinger equation 



d 



(5) 



(6) 



then follows from the variational principle SS — along 
with the boundary conditions <5|\l>(ti)) = 0; 5\^(t 2 )) = 0. 

Minimizing the action (|^) on a restricted variational 
basis for the wave function: 



dW* v V v = 1 



(7) 



leads to an effective action functional defined on the vari- 
ational parameters yi(t): 



dt ($ v \i^--H\iH v 



(8) 



with H given by 



(11) 



Given the above parametric form of the wave function, 
L(y,y) is always given by a function of the form, 



(12) 



where 



-t-oo 



i d 
,,!//) / dx-{V*(x;y) — V(x;y) 

d 



-V(x;y)—H>*(x;y)} (13) 



and 



%) 



dx^*(x;y)H^(x;y) 



(14) 



Minimization of the action, Eqn. (P-OP, leads to La- 
grange's equations: 



d_dL_ dL_ 

dt diji dyi 



= , for i = 1 , n . 



(15) 



The equations of motion for yi can be found easily using 
the specific Lagrangian defined in Eqn. (fl2|), 



Vj 



3=1 



dh{y) 
dyi 



(16) 



where the wave function is usually given in the coordinate 
representation. Extremization of the the effective action 
via 5r[?/i] = yields the dynamical equations obeyed by 
the variational parameters. 

In order to show that any variational solution leads to 
a symplectic Hamiltonian dynamics for the variational 
parameters (the case of Gaussians was considered in Ref. 
p3[), we consider general trial wave functions which are 
completely determined by n time-dependent functions of 
the form yi(t), i = 1, • • • , n, and written formally as 



V(x,t) =V(x; yi (t)) . 



(9) 



Here we choose for simplicity a one dimensional 
Schrodinger equation with arbitrary potential. Note that 
the entire time dependence of the wave functions is con- 
tained in the variational functions yi (t) . The Dirac form 
of the action is then given by 



r[y}= I dt I dx^*{x-y{t))\i--H\^{x-y{t)) 



d 



dtL(y,y) , 



(10) 



where My (y) is an anti-symmetric matrix given by 



Mij(y) 



dyj % 



(17) 



If the inverse of My exists, the equations of motion can 
be put in a symplectic form: 



N 



(18) 



Since My is also anti-symmetric, h(y) is a conserved 
quantity: 



dh(y) ^ dh . ^ dh x dh 



dt 



Following Das JTJ], we now introduce Poisson brackets 
by: 



I A B\-^ dA{y) M-' dB{y) 



(20) 
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So, for example, 



{y i ,y j } = Mr 1 



(21) 



The equations of motion can now be written in terms of 
these Poisson brackets: 



= {ftW=EV^ = DM}^ (22) 



the quantum mechanical system. At N = 1, the Gaus- 
sian approximation reduces to the well known Hartree 
approximation. 

The operator equations of motion following from the 
Hamiltonian (pq) are 



Xi + (m 2 + e 2 A 2 )xi = 
A + e 2 J2^A = 



(27) 
(28) 



The antisymmetry of the Poisson brackets is explicit in 
their definition (|2l|). However, they must also obey Ja- 
cobi's identity: 

{yi,{yj>yk}} + {y 1 Avk,y l }} + {y^iyuyj}} = o , (23) 

which is satisfied if My obeys Bianchi's identity: 

dM jk 



dM l3 dMki 



dyk 



dy dyi 



= 



(24) 



But Bianchi's identity is always satisfied for My of the 
form 



diiTj — dj%i . 



(25) 



Thus our definition of the Poisson brackets satisfies Ja- 
cobi's identity, and the set of classical equations of motion 
(p2|) are symplectic. 



III. HARTREE APPROXIMATION AND THE 
LARGE TV LIMIT 

We have shown that a time-dependent variational ap- 
proximation always leads to a Hamiltonian dynamical 
system for the variational parameters. Since such a sys- 
tem is generically nonlinear, there is a strong likelihood 
of chaos in the phase space of these Hamiltonian param- 
eters. In this section we derive two different approxima- 
tions for the coupled oscillator problem. The first keeps 
Gaussian correlations (Hartree approximation) for both 
oscillators, while the second (large N approximation) ig- 
nores fluctuations in the A oscillator. The second approx- 
imation has been derived previously from a path integral 
approach by making N copies of the x oscillator and 
then taking the large N limit. 

The model Hamiltonian that generalizes the two- 
oscillator problem to an N + 1 oscillator system is 



h = \a 



n N 

^ " ' .2 ' " , .. ■ 2 |2\ x 2 



^-pf + -(m' + e 2 A 2 ) 
i=i i=i 



(26) 



where we have introduced an N 



tor x^; ji = 0, 1 



1 component oscilla- 
, N with xa — A and the other N 
oscillators labeled by the roman indices i — 1, 2, ■ • • , N. 
We show below that at large N, a Gaussian ansatz for 
the wave function reproduces the exact large N limit of 



Taking expectation values of these two equations we ob- 
tain 



(xi)+m 2 {xi) +e 2 (A 2 x l ) = 
(A) +e 2 (x 2 A) =0. 



(29) 
(30) 



It was shown in Ref. Q that in the large N limit, fluctu- 
ations of the A oscillator are suppressed by 1/JV and the 
exact equations (p9|) and (J30|) are approximated by 



(xi) + m 2 (xi)+e 2 (A) 2 (xi) = 0, 
(A) +e 2 (x 2 )(A) =0. 



(31) 
(32) 



The semiclassical field (A) now has a time-dependent 
mass given by the quantum expectation value (x 2 ). The 
quantum oscillator Xi has a mass with time dependence 
controlled by (A). (This system was discussed in detail 
in Refs. It is also perfectly clear that the large N 

limit is equivalent to treating the A oscillator classically 
(i.e., ignoring the quantum fluctuations about the mean 
value of A). 

The equations governing (A) and (x 2 ) = G when 
(x) = were shown to be derivable Q from the effec- 
tive classical Hamiltonian: 

Heff = \v\ + 2^G + A + | ( m 2 + e 2 A2) q (33) 

We will show below that using a Gaussian trial wave 
function in Dirac's variational principle and taking the 
large N limit will lead to the same effective Hamiltonian 
( |33| ) for the evolution of the expectation values. However, 
if instead of taking the large N limit, we set N = 1, and 
a trial wave function which is a product of Gaussians in 
A and x, then the equations for the expectation values 
become: 



(x i )+m 2 (x i ) + e 2 (A 2 )(x i ) = 
(A)+e 2 (x 2 )(A) =0 



(34) 
(35) 



Here (A 2 ) = (A) 2 + D, and D is the Gaussian quantum 
fluctuation of the A oscillator (which also is the width 
of the A wave function) . In this case we will also get an 
effective Hamiltonian description of the dynamics, but 
with two more parameters, D and Hd- We will compare 
these two approximate Hamiltonian dynamics with the 
numerical simulation of the exact dynamics. 
Our choice for the trial wave function is 
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1 G 

V v (x M ) = Aexp[--(x - q(t))n(x - q(t)) v (— 01)^ 



+^P„(t)(x - q(t))„} 

where the normalization constant is given by 

1, 



(36) 



N 



exp 



-Trln(27r^G) 



The variational parameters are related to various expec- 
tation values taken with respect to the variational wave 
function 5'„: 



Pi{t) 

Gij(t) + qi{t)qj(t) 
2q i {t)p j (t) +4IL ik (t)G kj {t) 



(*„|Xi|*„) , 

(VvlxiPj+pjXilVv) . (37) 



The equations for these expectation values are obtained 
by varying the effective action, or equivalently from 
Hamilton's equations using the effective Hamiltonian. 

The effective action for the variational parameters 
p,q,G,U is 



r = I dt 



N 



J2PiQi +PaA - hTr[ilG] - H ef} 

,4=1 , 



(38) 



where TrfABl 



A^B,,^ and the effective Hamiltonian, 



H eff = (9 V \H\9 V ) 



N 2 2 
i=l 



2^Tr[nGn] 



iV 



[^- + ^(A 2 + G o)] 



i=l 



(qf + G u ) 



(39) 



This last equation gives the effective Hamiltonian for the 
dynamics of the N + 1 oscillators in the Hartree approx- 
imation. For simplicity (as was done in Ref. 0), we now 
specialize to the case q(t) = p(t) = 0. In this case G 
and n are diagonal (in general, they are also diagonal to 
leading order in the 1/N expansion). Since we have TV 
replicas of the x oscillator, the diagonal condition simply 
means that Gy(i) = G(t)Sij. Inserting this condition in 
(E39|) we find 



H 



(o) _ 



f 



eff-^PA + ^iNIl^G- 



U 2 D D) +-[— + - 



8 V G 



hN 



(A 2 + HD)] G. 



(40) 



Setting N = 1 in (ftO|), we find the effective Hamiltonian 
that controls the Hartree approximation: 



H 



(0) 

H 



1 



p 2 A + 2h(U 2 G G + U 2 D D) 



h r 2 

"2 [m 



(A 2 + hD)} G , 



(41) 



where G = Gu and D = Goo- 

Next we take the large N limit of (f40| ) using the same 
scaling argument as in determing the large N limit of 
the path integral formulation Qj: We let A — > \^NA 
and pa — > Pa (leaving invariant eA = eA). Dividing the 
effective Hamiltonian by AT and keeping the leading term, 
we find that the large N Gaussian effective Hamiltonian 
is exactly the same as the effective Hamiltonian found 
from the leading order large N action ||. The rescaled 
effective Hamiltonian reads 



H 



(0) _ h-(0) 
eff 



H eff/ N 



1 2 

= 2 Pa 



2hU 2 G G 



8G + 2^ 



e 2 A 2 ) G , (42) 



which is in complete agreement with ( J33| ) . (Tildes denot- 
ing the rescaled variables have been suppressed above.) 
At N — 1, the Hartree approximation has two more vari- 
ational parameters D and Hd compared to the large N 
approximation. These are related to the real and imag- 
inary part of the width of the wave function for the A 
oscillator and are obviously not incorporated in the large 
N approximation. Because of the extra degrees of free- 
dom incorporated in it, one might anticipate Hartree to 
be the better of the two approximations. 

In the Hartree approximation, the Hamilton's equa- 
tions for the expectation values are: 



A = pa , Pa = — e UAG , 

G = 4HU G G , D = 4HU D D , 



(43) 
(44) 



= ^2 ~ 2m l ~ \™ 2 ~ \^ ( A2 + hD ) > (45) 



n 



D = 



8G 2 
h 

8D 2 



- 2hn 2 n - -e 2 1l 2 G 



(46) 



In the leading order large N approximation, D — 0, and 
there is no equation for n^. 

For numerical work it is sometimes convenient to 
switch to a set of coordinates where the kinetic terms 
have the usual canonical form. Defining p 2 G — G and 
p D = D, the new Hamiltonian is 



im In 1 o 1 2 % f 1 1 

+ \[m 2 +e 2 (A 2 + h P 2 D )] p 2 G , 

with the resulting equations of motion 

A = pa , PA = -e 2 hAp G , 

Pg = Pg 1 Pd =Pd , 

p G = -^ T -(m 2 + e 2 h (A 2 + hp 2 D )) , 

_^ r 2fc2 2 

IPI 



Pd = — e 2 h 2 p D p G 



(47) 

(48) 
(49) 

(50) 
(51) 



Again the equations for leading order large N are ob- 
tained by setting pr> — and dropping po- The ad- 
vantage of this form is the ease in writing symplectic 
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integrators and also simplifying the form of the matrices 
needed to compute the Lyapvmov exponents. 

The above equations can now be solved numerically. 
Chaos (in the sense of nonzero Lyapunov exponents) ex- 
ists for large enough values of e 2 and for energy suffi- 
ciently above the ground state energy. 

IV. SEMIQUANTUM (GAUSSIAN) CHAOS 
A. Numerical Methods 

In this section we display evidence that both the large 
N and Hartree approximations are chaotic for appro- 
priate values of the energy E and the coupling e. (In 
Ref. Q, the large N approximation alone was shown to 
be chaotic.) The dynamics of test trajectories in the 
above approximations was studied using a fourth-order 
symplectic integrator. (This integrator was implemented 
using the second set of variables defined at the end of 
the last section.) Chaos was characterized quantitatively 
by measuring the Lyapunov exponent for different initial 
conditions using standard techniques [jl5| . 

In order to check whether the chaos seen in the ap- 
proximation is of some relevance to the full quantum 
problem, a numerical solution of the corresponding time- 
dependent Schrodinger equation is required. This was 
accomplished by using second and fourth-order unitary, 
split-operator, spectral solvers that we have recently im- 
plemented on a large parallel computer p6| . By us- 
ing large grids (up to 4096 x 4096) sufficient resolution 
is achieved to accurately evolve the wave function over 
times long enough to meaningfully compare with results 
from the variational approximations. 

The phase space of the large N and Hartree approx- 
imations was characterized using Poincare sections. At 
relatively low energies and modest values of the coupling 
constant e, both the approximations led to integrable dy- 
namics. However, increasing either the energy or the cou- 
pling constant quickly led to nonintegrability. While not 
carrying out an exhaustive analysis, we did isolate para- 
metric regions where the chaos was relatively soft (the 
area of stochastic orbits was small compared to the area 
occupied by regular orbits) and regions where the dynam- 
ics was predominantly chaotic. We also ran a large set of 
initial conditions to sample the regions in coupling con- 
stant/energy space where the approximations were reg- 
ular. This was accomplished by implementing a parallel 
code to compute the Lyapunov exponents for a large set 
of independent trajectories. 

B. Numerical Results 

There are two separate but related questions concern- 
ing the variational approximations. The first question re- 
lates to how well they track the exact numerical calcula- 



tions. We find that the approximations break away from 
the exact calculation on a short time scale independent 
of whether they are chaotic or not. However, qualitative 
agreement with the numerical results is much better in 
the nonchaotic case. The second question refers to the 
stability of the approximate solutions as well as the ex- 
act solution. In the chaotic regime of the approximations, 
the approximate evolution is sensitively dependent on ini- 
tial conditions whereas the exact evolution is not. After 
a finite time, two approximate evolutions starting from 
almost identical initial conditions become completely dif- 
ferent in the chaotic case and no longer bear any phase 
relationship amongst themselves or to the exact solution. 
This is in contrast with the behavior in the integrable 
case. 

The addition of variational parameters has two effects: 
it qualitatively improves the long time behavior in both 
the regular and chaotic regimes even though the break 
time from the exact behavior is not affected. Secondly, 
there is some evidence that the onset of chaos is de- 
layed as more parameters are added and that the value 
of the maximum Lyapunov exponent is also decreased. 
However, an exhaustive study would require a system- 
atic method of adding variational parameters for the trial 
wave functions and this we leave to the future. 
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FIG. 1 Domain of integrability for the large N approxi- 
mation in energy and coupling constant space. The phase 
space was sampled with ten initial conditions at each (e, E) 
point (and trajectories asymptoting to positive Lyapunov ex- 
ponents were searched for. At fixed e, the region above any 
point denoted by the top square in the figure corresponds to 
chaotic dynamics, i.e. in the set of trajectories sampled there 
was at least one with asymptotically positive Lyapunov ex- 
ponent. The region below the bottom square corresponds to 
integrable dynamics. 

In Fig. 1 the approximate region of regularity for the 
large iV approximation is displayed. Each (e, E) point 
was sampled by ten trajectories, with the Lyapunov ex- 
ponent calculated for each. Within the uncertainties of 
our sampling scheme the integrable and nonintegrable 
regions cannot be sharply distinguished: the top set of 
points denotes at least one trajectory having an asymp- 
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totically positive Lyapunov exponent while below the 
bottom set of points no such trajectory was ever found. 
The true boundary should be roughly in the middle of 
these two curves. The results for the Hartree approxima- 
tion are very similar and slightly above the integrability 
curve for large N but the difference is of order the un- 
certainty band. Whether there is a general (monotonic) 
tendency for this to happen as the number of degrees of 
freedom is further increased is an interesting speculation 
which needs to be explored further. 
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FIG. 2 A typical computation of the maximal Lyapunov 
exponents for the large N (upper curve) and Hartree (lower 
curve) approximations. Parameter values for this run were 
e = 1 and E = 5. 

The Lyapunov exponents for the two approximations 
were computed in the chaotic parameter regime. For all 
cases we studied the maximal exponent in the Hartree 
approximation was less than the corresponding exponent 
in the large N approximation. A typical example of these 
results is given in Fig. 2. 

Poincare sections are another way to explore the do- 
mains of integrability for the two approximations. For 
the "boundary" regions of Fig. 1, the phase space was 
largely mixed, with stochastic regions coexisting with 
regular regions. We checked for random values of the 
parameters that the region below this boundary was reg- 
ular. Above, it was dominantly chaotic. It was difficult 
to use Poincare sections for the Hartree approximation 
because more degrees of freedom means running much 
longer to get acceptable statistics. We did run checks for 
a few parameter values and found results consistent with 
Fig. 1 including the fact that chaos occurred at larger 
values of the parameters. For example, while the large N 
approximation had appreciably chaotic regions at E = .8, 
e = .7 the Hartree approximation was completely intc- 
grable for those values of the parameters. (Note that the 
energy E is different for the large N and Hartree ap- 
proximations since D and Hu contribute in the Hartree 
approximation, but not in large N .) For the parameter 
values, e = 1 and E — .8 we show two Poincare sections 
in Figs. 3 and 4 (large N) which are typical for values of 
the parameters near the boundaries of Fig. 1. 
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FIG. 3 A Poincare section in the A, A plane for e = 1 and 
E — 0.8 for the large N approximation. The phase space was 
sampled by 256 different trajectories. 
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FIG. 4 A Poincare section in the G, Ug plane in the large 
N approximation for the same set of parameters as Fig. 3. 

In order to assess the relevance of the chaos seen in the 
approximations we have compared the approximate evo- 
lutions with exact numerical solutions of the Schrodingcr 
equation with Gaussian initial data. The exact evolu- 
tion shows no hint of the sensitivity to initial conditions 
exhibited by the approximate dynamics. As illustrated 
in Figs. 5, 6, and 7, in both the regular and chaotic 
regimes the approximations quickly deviate from the ex- 
act results on a time scale of order unity, this signaling 
the breakdown of the Gaussian approximation. 

The Lyapunov time sets a maximum time for which 
the approximations can agree with the exact quantum 
dynamics. In fact, consistent with this statement we ob- 
serve that the time of breakdown of the approximations 
and the Lyapunov time are of the same order. However, 
this should not lead one to conclude that the accuracy 
dramatically improves when the approximate dynamics is 
integrable. Indeed, even in integrable parameter regimes, 
the breakdown time can remain of order unity (Fig. 5). 
Therefore, for coupling constants of order unity, these ap- 
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proximations tend to be rather poor. This is because sig- 
nificant non-Gaussian structure forms in the exact wave 
functions relatively rapidly. 
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FIG. 5 Evolution of (A) for e = .3 and E = 1. This is 
within the parameter range for nonchaotic evolution within 
the approximations. Both approximations break away from 
the exact evolution at t ~ 5 but stay in phase at later times. 
The Hartree approximation (H) does better in tracking both 
phase and amplitude. The curve with the smallest average 
amplitude corresponds to the exact quantum evolution (Q), 
and the one with the largest average amplitude corresponds 
to the large N expansion (N). 
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FIG. 6 Evolution of (A) for e = 1 and E — 5. The approx- 
imate evolutions are now chaotic. They break away from the 
quantum evolution at time t ~ 2 (denoted by the point 1 in 
the figure) and break away from each other at point 2 (t ~ 4) . 
In this case the evolutions quickly dephase from each other 
and from the quantum evolution. 

The chaos inherent in the approximations is demon- 
strated in Figs. 8, 9, and 10, for the evolution of (A), G, 
and D. In these figures we show two trajectories for each 
of the approximations, one corresponding to an initial 
G = .5 and the other to G = .5001 (all other parameters 
held fixed). The deviations of these two curves are con- 
sistent with the calculated Lyapunov exponent (which is 
of order unity) and an initial deviation of order 10 A 















A ; 






Large N / 




/ \ 
/ Hartree,-.! 




1 


-?-—// 


V.. -•■--/..Quantum \ 





2 4 6 8 10 



FIG. 7 Evolution of G for the same parameters as Fig. 6. 
The break from the quantum evolution occurs at t ~ 1 (de- 
noted by the point 1 in the figure) and the approximations 
break away from each other at point 2 [t ~ 2). Even though 
the Hartree approximation is not correct it does not have the 
big excursions shown by the Large N approximation. 
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FIG. 8 Evolution of (A) for the same parameters as Fig. 
6. Points 1 and 2 mark the breaking away of two nearby tra- 
jectories in the Hartree and Large N approximations. After 
this time, the trajectories rapidly dephase from each other. 
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FIG. 9 Evolution of G for the same parameters as Fig. 
6. The trajectory denoted by Q is the quantum evolution. 
Points 1 and 2 mark the breaking away of two nearby trajec- 
tories in the Hartree and Large N approximations. 
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FIG. 10 Evolution of D for the same parameters as Fig. 
6. The quantum evolution and the Hartree approximation 
deviate from each other at t ~ 1 and the two nearby trajec- 
tories of the Hartree approximation break from each other at 
t ~ 15. 

The approximations discussed here break down when- 
ever there is significant non-Guassian structure in the 
actual wave function. As long as the coupling is of or- 
der unity this happens relatively rapidly Examples of 
the numerically evaluated probability densities in A are 
shown in Figs. 11 and 12 for for values of the param- 
eters which correspond to integrable and nonintegrable 
evolutions. 
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FIG. 11 The initial and final (t = 100) probability densi- 
ties for A with e = .3 and E — 5 (integrable case). 
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FIG. 12 Same as Fig. 10 with e = 1, E = 5 (nonintegrable 
case). In this case the final time, t — 40. 



V. CONCLUSIONS 

The central results of this investigation may be en- 
capsulated succinctly: all time-dependent variational ap- 
proximations based on a Dirac approach are Hamiltonian 
and generically nonlinear. Therefore all such approxi- 
mations can be chaotic. Since exponential divergence 
of expectation values in time is ruled out in full quan- 
tum mechanics, the Lyapunov time associated with the 
approximate evolution sets a time scale beyond which 
the approximation breaks down. We have investigated 
this last point in two particular examples (large N and 
Hartree for a two-dimensional potential), find both to be 
chaotic, and by comparison against numerically obtained 
solutions, show explicitly that the approximations break 
down on the Lyapunov time scale. We also show that 
even in nonchaotic regimes, the approximations break 
down very quickly. Thus the mere absence of chaos is not 
an indicator of the accuracy of these approximations. 

We would also like to point out that suggestions have 
been made in the literature that semiquantum chaos may 
in fact be a real effect (e.g., Ref. H and rather more 
strongly in Ref. fl7|| ). However, these claims were not 
backed up by careful comparisons with exact calcula- 
tions. The detailed results reported here, along with 
the fact that Gaussian approximations are dynamically 
completely classical ||, imply exactly the opposite con- 
clusion (in substantial agreement with the arguments of 
Ref. @). 

The fact that in the chaotic regime, the approxima- 
tion signals its own breakdown has an interesting phys- 
ical consequence: if the 1/N approximation is in fact 
sensible then a breakdown at leading order must imply 
that the next-to-leading terms are becoming large on the 
same time scale. Since, in field theory the leading order 
approximation does not incorporate collisions, what this 
implies is that the collisional time scale can be estimated 
from the breakdown of the leading order result itself, 
without actually having to compute the next-to-leading 
order contribution. Given the complexity of higher-order 
calculations this feature may be extremely useful. This, 
and other aspects of the field theoretic problem are now 
under investigation. 

One way to of incorporating higher order correlation 
functions in dynamical approximations is to consider trial 
wave functions of the form Gaussian times polynomials. 
This can be put in correspondence with the large N ex- 
pansion which can be shown to lead to the same struc- 
ture. An interesting question is whether opening up the 
possibility of including higher order correlations in this 
way will improve the long time behavior of the variational 
approach. 
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